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ABSTRACT 

The Cosmic Background Explorer (COBE) spacecraft will exhibit complex 
attitude motion consisting of a spin rate of approximately -0.8 revolution 
per minute (rpm) about the x-axis and simultaneous precession of the 
spin axis at a rate of one revolution per orbit (rpo) about the nearly per- 
pendicular spacecraft-to-Sun vector. The effect of the combined spinning 
and precession is to make accurate attitude propagation difficult and 
the 1-degree (3 a) solution accuracy goal problematic. 

To improve this situation, an intermediate reference frame is introduced, 
and the angular velocity divided into two parts. The nonspinning part is 
that Avhich would be observed if there were no rotation about the x-axis. 

The “spinning” part is simply the x-axis component of the angular veloc- 
ity. The two are propagated independently and combined whenever the 
complete attitude is needed. This approach is better than the usual one- 
step” method because each of the two angular velocities look nearly con- 
stant in their respective reference frames. Since the angular velocities are 
almost constant, the approximations made in discrete time propagation 
are more nearly true. 

To demonstrate the advantages of this “nonspinning” method, attitude is 
propagated as outlined above and is then compared with the results of the 
one-step method. Over the 100-minute COBE orbit, the one-step error 
grows to several degrees while the nonspinning error remains negligible. 


"This work was supported by the National Aeronautics and Space Administration (NASA) /Goddard 
Space Flight Center (GSFC) , Greenbelt, Maryland, under Contract NAS 5-31500. 


PRECEDING PAGE BLANK NOT FILMED 



COBE ATTITUDE AND THE NEED FOR MORE ACCURATE 

PROPAGATION 


The attitude of the Cosmic Background Explorer (COBE) is three-axis stabilized with the 
minus x-axis maintained 94 degrees from the Sun line (Figure 1). The spacecraft pitches 
about that line at a rate of 1 revolution per 100-minute orbit period. It also spins about 
the minus x-axis once every 75 seconds. COBE attitude can be expressed as an 
Euler 3-2-1 (pitch-roll-yaw) rotation sequence with respect to a rotating Earth-Sun coordi- 
nate frame with the z-axis pointing toward the Sun and the y-axis pointing along the cross 
product of the Sun and Earth vectors. Nominal pitch and roll are then 0 and -4 degrees, 
respectively. 



Figure 1. COBE Attitude Profile 


Propagating the attitude from samples of the angular velocity assumes that the angular 
velocity remains constant in the body frame over the sample interval. The COBE gyro 
sampling intervals are 0.5, 1, 2, or 4 seconds, depending on the telemetry format and 
data rate. For COBE, which spins as much as 19.2 degrees per gyro sample, the pitch 
component of the angular velocity can change direction significantly. Under these circum- 
stances, the usual method of propagation, which handles the total incremental rotation at 
once, introduces errors that accumulate over time and become unacceptably large. 

This paper describes a variation on the usual one-step propagation that reduces this error 
by introducing an intermediate nonspinning coordinate frame. The advantage of the 
nonspinning frame is that the angular velocities used for its propagation and for the 
subsequent transformation to the body frame vary much less between gyro samples. The 
associated equations for gyro calibration follow along with numerical estimates of the 
improvement in COBE propagation accuracy. 
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NOTATION AND THE USUAL ONE-STEP ATTITUDE PROPAGATION 

Attitude is represented here by the orthogonal inertial-to-body coordinate transformation 
matrix, A B /i , and the kinematic equation for its propagation (Reference 1, p. 512) is 

Ab/i = ~ A B /i (1) 

The dot (•) above Ab/i indicates differentiation with respect to time. Since the angular 
velocity, (Db , is not constant, Equation (1) is solved numerically. Still, a formal solution 
may be written as 


Ab/i (t) = ^b/i (t» to) Ab/i (to) 

<J>b/i is the attitude propagation matrix satisfying the differential equation: 

<I>B/I = ~ {*?B ^B/I 

and the initial condition: 


^b/i (to, t 0 ) = I 


(2) 

(3) 

(4) 


where I is the identity matrix. 

The subscripts I and B refer here to the inertial and body coordinates in which a vector or 
matrix is expressed. The vector o> B , for example, is the angular velocity in body coordi- 
nates, and (o B is the antisymmetric matrix derived from it. 
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<^B1 
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(5) 


The slash (/) indicates a transformation from the frame on the right to that on the left. 
Thus, A B /i is the inertial-to-body coordinate transformation matrix. 

NONSPINNING INTERMEDIATE FRAME AND TWO-STEP PROPAGATION 


In body coordinates, the spin component of the COBE angular velocity, <o'b , is constant, 
while the pitch component, co'b > varies with time: 


(Ob (t) = (o'b (t) + (o'b (6) 

If the spin axis is denoted by §b> these two portions of the angular velocity can be com- 
puted as follows: 

co'b - (I - ^ b (s b ) T )a>b (2) 
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COfi - S B ( S b ) T 0) b (8) 

An intermediate, nonspinning coordinate system denoted by the subscript N can be intro- 
duced that is defined by the propagation equation 

An/i = ~ A N /i (9) 

where 

o»n = (Ab/n) T (10) 

Since the magnitude of (o' B is less than that of a) B , a>' B does not change as much as (o B 
does between gyro samples. Thus, all other things being the same, the error in propagat- 
ing A N /i should be less than the error in propagating A B /i . If the propagation from the 
nonspinning frame to the body frame can be done perfectly, as is plausible since the spin 
rate is constant, the total propagation error for this two-step method should also be less 
than that for the usual one-step method. 

To complete the propagation of A b/i , it remains to compute A B/N , which transforms 
from the nonspinning to the body frame: 


A B /n = A B /[ (An/i) t (11) 

By the product rule for differentiation, the A B/] is 

Ab/n = A B /i (A n/ i) t + A B /j (An/i) t (12) 

Substituting for A B /i and A N /i from Equations (1), (6), and (9) and combining attitude 
transformations yields 

A B /n = - (^b + (|)b) A B /n + A B /n &n (13) 

Noting that is the similarity transformation of w B , 

= (A B /n) t w b A b/n (14) 

gives the differential equation for the propagation of A B / N : 

A B /n = - A B /n (15) 

With equations for both A b/n and A n/] , the complete attitude, A b/i , can be propagated. 

GYRO CALIBRATION 

In addition to reducing the propagation error due to the finite gyro sampling interval, it is 
usually necessary to calibrate the gyros to reduce systematic errors in the sensed angular 
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velocity. Because calibration involves the same type of computations as propagation, it is 
also done more accurately with a two-step method. 

With one-step propagation, gyro calibration errors, A a , can be found using the solution 
to the following error propagation equation (Reference 2, p. 4-13): 

9 B = -w B 9 B + Aw B (16) 

where 9 B is the attitude error expressed in axis and angle form, and the angular velocity 
error, A a> B , is related to the gyro calibration errors by the matrix G(<u B ). 


Aw B = G((»b) Aa (17) 

The solution to the error equation has the form 


e B (t) = r B/ , (t, t 0 ) Aa + 0> B/ I (t, to) 0 B (to) (18) 

where the variational matrix, r B /i , transforms the gyro calibration errors into contribu- 
tions to the attitude error 

r B/I (t, t 0 ) = f <D b/1 (t,r) G (<u B ) dr (19) 

■'to 


The matrix r B/ i then serves as the partial derivative of the propagated attitude error with 
respect to the gyro parameter errors: 


d9 B 
dAa 


= Tb/I (t, to) 


( 20 ) 


If one knows 0 b at times t and t 0 , A a can be found from Equation (18). 

A corresponding variational matrix is needed for calibration with the nonspinning propa- 
gation method. Because the two methods are different, there is no reason to expect the 
variational matrices to be the same. Both steps of the nonspinning propagation, however, 
follow the same kind of propagation equation as does the one-step method, and the corre- 
sponding error equations can be applied to each step separately: 


= Ao>n - 0n (21) 

0 B = Aft> B - o) B 0 B (22) 

Here, 9'u is the error in the nonspinning propagation expressed in the nonspinning frame. 
9 b is the corresponding error in the spinning propagation expressed in body coordinates. 
The angular velocity errors, Acon and Acd b , are defined as follows: 

A^n = A n/b [I - §B (s b ) T ] G (g>b) Aa (23) 
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Ao>b = (§b) t G (<y B ) Aa (24) 

and O'b can be solved for as in the one-step propagation to give solutions of the form 

= r N /i A a + On (t c ) 

= r B /N A a + 0 B /N #B (to) 

where the variational matrices Tn/i and T B /n are computed as follows: 


Tn/i = J ^n/i (t, r) A<Wn dr 

■'to 

(27) 

Tb/n = f ^b/n (t> r) A<y B dr 
■'*0 

(28) 

The total propagation error, Os, then equals 


$b = A B /n #n + O'b 

(29) 


This gives the partial derivative of the attitude error with respect to the gyro calibration 
errors as 


(25) 

(26) 


Tb/I - r B /N + ^B/N Tn/I (30) 

which can be used in Equation (18) to solve for Aa. 

NUMERICAL SIMULATION 

Although the large COBE gyro sampling interval can be expected to degrade one-step 
propagation accuracy, it is useful to know how much of an effect it actually has. The 
nonspinning method must justify its additional computation with significantly better accu- 
racy. 

To compute the propagation error for each method, the pitch rate, 0, roll rate, <9, and 
spin rate, 0 , measured with respect to the inertial frame are assumed to be constant. 


0 = -360°/6000 sec = -0.06°/ sec 

(31) 

0=0 

(32) 

0 = - 0.8 rpm = -4.8°/ sec 

(33) 
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The exact attitude is then found by first performing the pitch rotation, 0 , followed by the 
roll and spin rotations, 6 and 0 . 


0 = 

(t - to) 0 + 00 

(34) 


o 

1 

II 

(35) 

0 = 

(t - to) Ip + 00 

(36) 


The one-step attitude can be propagated from the following formula for the angular veloc- 
ity (Reference 1, p. 765): 

cos (t) T = [-ip - sin (0) 0, cos (0) sin (0) 0, cos ( 6 ) cos (0) 0] (37) 

The two angular velocities for the nonspinning propagation are as follows: 

(y^ (t) T = [0, cos (cueit) coq 2 ~ sin (wbiO a>B3, sin (ft>Bit) <Wb 2 + cos (<WbiO &>b3] (38) 

o>b T = [0 - sin ( 6 ) ip, 0, 0] (39) 

The one-step and nonspinning attitudes are then compared to the exact attitude, and the 
angular differences are computed. A plot of the one-step propagation error for 0.5-second 
sampling intervals is shown in Figure 2 for a timespan of one orbit. While the one-step 
error grows to a maximum value of 2.4 degrees, the nonspinning error remains less than 
0.003 degree. 

The results show that, even for the smallest gyro sampling interval, which is 0.5 second, 
the one-step propagation errors are quite large. This is counterintuitive. Since a constant 
spin can be propagated without error, it would be expected that adding a much smaller 
constant pitch rate would have a negligible effect. The flaw in that argument is that, 
although the spin is constant in the body frame, the pitch angular velocity is not. It is 
constant in the inertial frame. The assumption about the size of the error is also slightly 
misplaced. Small pitch rates do produce slowly growing propagation errors, but because 
the orbit is correspondingly longer, they have more time to grow. 

Even more surprising than the size of the error is its oscillation. Rather than grow without 
limit, the propagation error peaks at the middle of the orbit. The reason is that while the 
actual pitch angular velocity moves continuously in the body frame, the sampled pitch 
angular velocity is fixed in the direction it has at the start of the interval. Thus, the 
sampled value lags the true value and introduces a roll component of angular velocity 
(Figure 3). Over an orbit period, the roll direction changes by 360 degrees, and the 
propagation error that builds up in the first half of the orbit decreases over the second 

half. 
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Figure 2. One-Step Propagation Error Time Dependence 



Figure 3. Sampling introduces a Roll Angular Velocity 


SUMMARY 


The nonspinning propagation method described here is a means of trading computation 
for accuracy when the body angular velocity changes direction between gyro samples. 
This method is currently implemented for batch attitude determination in the COBE flight 
dynamics support system (Reference 3, p. 3.1.2.42-1), where it is needed to meet the 
attitude determination accuracy goal of 1 degree (3o). Further investigation should still 
be done on understanding the effects of nonconstant roll, interpolating the angular 
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velocity to the midpoint of the sampling interval, and using higher order numerical inte- 
gration methods. 

Whether nonspinning propagation is worth the extra work for other missions depends on 
the magnitude and form of the angular velocity, the gyro sampling rate, and the accuracy 
requirements. The unexpectedly large errors that would have been observed for COBE, 
however, argue for consideration of this effect whenever a three-axis stabilized spacecraft 
undergoes a compound rotational motion, such as spinning and pitching. 
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